Crossover from Luttinger liquid to Coulomb blockade regime in carbon nanotubes 
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We develop a theoretical approach to the low-energy properties of ID electron systems aimed 
to encompass the mixed features of Luttinger liquid and Coulomb blockade behavior observed in 
the crossover between the two regimes. For this aim we extend the Luttinger liquid description 
by incorporating the effects of a discrete single-particle spectrum. The intermediate regime is 
characterized by a power-law behavior of the conductance, but with an exponent oscillating with the 
gate voltage, in agreement with recent experimental observations. Our construction also accounts 
naturally for the existence of a crossover in the zero-bias conductance, mediating between two 
temperature ranges where the power-law behavior is preserved but with different exponent. 
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The recent progress in nanotechnology has allowed 
for a detailed study of the transport properties of one- 
dimensional (ID) electron systems. The discovery of 
carbon nanotubes (CNs) in 19910, as a by-product of 
fuUerene production, has opened a new field of research in 
mesoscopic physics, especially because of their potential 
application to nanoelectronic devices 0- 

CN are rolled up sheets of graphite forming tubes that 
are only nanometers in diameter, with metallic or semi- 
conducting properties depending on the roll-up direction 
of the sheet. Regarding the metallic CNs, their transport 
properties are dictated by the low dimensionality of the 
system, which gives rise to strong electronic correlations. 
Thus, the effects of the Coulomb interaction become very 
important in CNs, and the way they manifest themselves 
depends mainly on the size of the system, the tempera- 
ture, and the quality of the contacts used in experiments. 

When the contacts are not highly transparent, the 
measurements of the conductance and the differential 
conductivity reflect the strong Coulomb repulsion in 
CNs. For temperatures that are typically below 1 K, 
the zero-bias conductance shows oscillations character- 
istic of the so-called Coulomb blockade regimeQ- The 
gate voltage between two peaks is related to the energy 
required to overcome the Coulomb repulsion when adding 
an electron between the barriers created by the contacts. 
In general, the intensity plots in terms of the bias and 
gate voltages show in that regime a sequence of diamonds 
where the electron tunneling is suppressed, with the size 
of the diamonds giving a measure of the energy needed 
to add an electron to the system^. 

For higher values of the temperature, such that the 
thermal energy is much larger than the level spacing, 
the transport measurements reflect instead the many- 
body properties of the system. The conductance shows 
then a different kind of suppression, which arises from 
the strongly correlated character of the electron liquid 13]. 
In a ID system, the electron quasiparticles are unstable 
under the slightest interaction and the low-energy excita- 
tions take the form of charge or spin density fluctuations. 



In the so-called Luttinger liquid regime, the absence of 
electron quasiparticles gives rise to a power-law suppres- 
sion of the tunneling density of states (TDOS) near the 
Fermi level, n(e) e", which translates in turn into a 
similar behavior of quantities like the zero-bias conduc- 
tance as a function of the temperature, G{T) ~ T". 

Such kind of signatures of Luttinger liquid behavior 
have been observed in single- walled as well as in multi- 
walled nanotubes (MWNTs)0- It has been shown that 
the zero-bias conductance and the differential conduc- 
tivity follow closely a power-law behavior as functions of 
the temperature and the bias voltage, respectively. In the 
case of the zero-bias conductance, such power-law depen- 
dence holds over a wide range of temperatures, down to 
values where the discrete character of the spectrum starts 
to become manifest. It is observed then the transition to 
the Coulomb blockade regime, where the suppression of 
the conductance arises from the single-electron tunneling 
mentioned above. 

At this point, it would be desirable to have a theoret- 
ical description of the crossover as the temperature de- 
creases from the Luttinger liquid regime to the Coulomb 
blockade regime. _There have been experiments like that 
reported in Ref. ^ in which the intermediate regime has 
been explored, measuring the zero-bias conductance at 
temperatures where the thermal energy becomes compa- 
rable to the level spacing in the discrete spectrum. The 
results reported in that reference provide evidence for the 
power-law behavior of the TDOS in MWNTs, but with 
an unusual dependence of the exponent a on the gate 
voltage. Clear oscillations have been observed in that 
plot, with a characteristic scale of the order of 1 V. An- 
other remarkable observation has been the existence of 
an inflection point in the log-log plot of the conductance 
versus temperature, for gate voltages corresponding to 
peak values in the a exponent^. 

The aim of the present paper is to provide a frame- 
work that may encompass the mixed features of Luttinger 
liquid and Coulomb blockade behavior observed in the 
crossover between the two regimes. For that purpose we 
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will rely on the ID many-body approach'^ , incorporating 
at the same time the discrete character of the spectrum 
which arises from the finite size of the system. Thus we 
will see that the polarizability, which is proportional to 
the single-particle density of states, is modulated at suf- 
ficiently low temperature due to the peaks building the 
density of states in the discrete spectrum. Such a modu- 
lation is translated to the scaling behavior of the quasi- 
particle weight, giving rise to oscillations in the exponent 
of the TDOS as the Fermi level of the system is shifted. 
We will show that the period estimated for the oscilla- 
tions is in agreement with the experimental observations, 
and that the proposed framework is able to account for 
the features seen in the measurements of the zero-bias 
conductance reported in Ref. '3. 

In our approach, a central role in distinguishing the 
Coulomb blockade from the Luttinger liquid regime is 
played by the energy distribution function nle^ep)- At 
low temperatures, that function has peaks that reflect 
the level spacing Ae arising from the level quantization 
in a finite nanotube. We can model each peak with a 
normalized function /(e,n), with a shape depending on 
the temperature T 



fcsTcosh^(^) 



(1) 



The energy distribution function is cutoff by the Fermi 
energy ep, as it is given by 



ep/Ae 

i(e;e_F) = ^ /(e. 



(2) 



n=0 



At sufficiently high temperatures, we can neglect the 
level spacing and calculate the energy distribution func- 
tion by integrating 



n(e;e_F) 



eF/Ae-l/2 



dyf{e, v) 



1 + g(£-eF + Ae/2)/fcBT ' 



(3) 



The last formula gives the usual Fermi-Dirac distribution. 
We realize however that the energy distribution function 
has in general a more structured shape, that depends on 
the ratio between the level spacing Ae and the thermal 
energy ksT as illustrated in Fig. ^ 

The structure of peaks in the energy distribution func- 
tion affects several objects in the many-body theory. This 
is the case of the one-loop polarization n(°)(g,a;q), which 
is given in terms of the electron propagator G^^\k,uJk) 
by 



(4) 

This object counts the number of particle-hole excita- 
tions that can be built with momentum q across the 




FIG. 1: Plot of the energy distribution function n{e\eF) 
for different values of the temperature corresponding to 
ksT / Ae — 0.25,0.6, 1.0. The solid lines have been obtained 
from Eq. (|5J and the dashed lines by using the Fermi-Dirac 
distribution @. 



Fermi level. Thus, in the case of a model with linear 
dispersion e{k) « vpk, it becomes proportional at zero 
temperature to 



dk 0{kF - k)9{q + k - kp) = q. 



(5) 



However, in the model with a distribution function given 
by (O, the peaks give rise to a factor in the integration 
that reflects the discrete structure 



dk 'n(vFk]eF) — n{vp{q + k)]eF)) — qNT{sF)- 

(6) 

The usual form of the ID polarization is modified then 
by the presence of the last factor in © : 
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(7) 



with the factor Ns standing for the number of subbands 
crossing the Fermi level. The factor Nt{sf) is itself a 
periodic function of ep, since the position of the Fermi 
energy fixes the portion of the last peak that has to be 
integrated in the spectrum, in order to obtain the polar- 
izability. 

We can incorporate the polarization in the RPA 
evaluation of the electron self-energy I](fc, w^). This can 
be expressed as 



Y.{k,u!k) = i 



A/VF 



dq 
2^ 



V{q) 



dujq 

"2^ 



l-V{q)Iim{q,LJg)-' 



(8) 



where V{q) stands for the ID Coulomb potential in mo- 
mentum space. The high-energy cutoff A comes from the 
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existence of a microscopic short-distance scale in the elec- 
tron system and, for our purposes, it provides a way of 
having under control the scale dependence of the different 
observables0. 

The RPA scheme allows for an accurate description of 
the behavior of the quasiparticle weight at low energies, 
since it gives the exact result for the anomalous dimen- 
sion of the electron field in the ID system 10]. Thus, the 
self-energy turns out to have terms linear in uik that de- 
pend logarithmically on the cutoff A. This is the signal 
of the energy dependence of the scale Z^/^ of the electron 
wavefunction. The dressed electron propagator G{k,ujk) 
is given by 

+Z-\K){ujk - VF(Txk)-i{g)\og{K) (9) 

where g is the effective coupling constant, given in terms 
of a suitable average value of the Coulomb potential V 
hy g = 2NsNT{eF)V /nvp. The function j{g) bears also 
an explicit dependence on the oscillating factor iVj-(e^), 
as it is given by 

The requirement of cutoff-independence of the renor- 
malized Green function G is what forces the dependence 
of Z on the energy scale A, as it is bound then to satisfy 
the flow equation 

A-^logZ{A)=j{g). (11) 

The quantity at the right-hand-side of (|ll|l is actually the 
exponent a of the power-law behavior characterizing the 
suppression of the quasiparticle weight in the low-energy 
limit A — > 0. The novelty of our framework is that such 
an exponent depends now on the Fermi energy of the 
system, as the factor A't(£f) is sensitive to the position 
of the Fermi level in the structure of single-particle levels. 

We now ascertain that our approach may account for 
the experimental features reported in Ref. @. For this 
purpose, we have first to check that the Fermi energy 
shift As producing an oscillation in the a exponent is 
consistent with the variation of the gate voltage applied 
in the experiment. The dependence of a on the gate 
voltage Vg has been plotted in Figs. 2(c) and 2(d) of 
the mentioned reference. There it can be seen that, in a 
range of 20 V for the gate voltage, there are about eight 
oscillations in the a exponent. The corresponding change 
in the Fermi energy can be inferred from the capacitive 
coupling between the gate and the MWNT, that has been 
estimated to be about 1 x 10"^ • Then, each oscillation 
in the plot should correspond to a variation in the Fermi 
energy of about 2 meV. This matches well the value of Ae 
estimated by taking the level spacing from the finite nan- 
otube length L as a lower bound, Ae ~ hvp/L ~ 1 meV, 
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FIG. 2: Plot of the exponent a as a function of the Fermi en- 
ergy ef (in units of Ae) for different values of the temperature 
T= 1,40, 150 K. 

hence providing a strong indication that the quantization 
of the nanotube levels is at the origin of the oscillations 
observed in the a exponent. 

The power-law behavior obtained from Hll|l applies di- 
rectly to the TDOS n{e), since we have 

n{e) NsZ{e) ^ e^^^l (12) 

We can compare the dependence of the exponent 7(3) 
on Ep with that of the power-law behavior of the con- 
ductance measured in Ref. ^ In our framework, the 
exponent depends on the position of the Fermi level due 
to the modulation introduced by Nt{£f) in 7(5)7 but also 
because of the possible variation of the Fermi velocity vp 
as the Fermi level is shifted. In this respect, we recall 
that MWNTs are usually hole -dop ed, with several sub- 
bands crossing the Fermi level[l3, so that a significant 
decrease in the Fermi velocity at some of the Fermi points 
is likely as the top of the lowest subband is approached. 
Actually, a drift in the average value of the exponent a 
is appreciable in Figs. 2(c) and 2(d) of Ref. 'S, as the 
gate voltage runs from —10 V to 10 V. We have incor- 
porated this effect by taking into account a suitable shift 
of the Fermi velocity in our model, from vp to hvp/Q, as 
£p runs over 10 oscillations in the exponent. Finally, we 
have taken for the interaction strength V ~ bvp, bear- 
ing in mind the screening from the metallic shells of the 
MWNT 131]. The corresponding plot for a as a function 
oi £p is given in Fig. |21 We observe the similarity of 
the plot with the experimental curves of Ref. i8j for the 
exponent over the 20 V range in the gate voltage. 

Our approach provides also information about the de- 
pendence of the exponent a on the temperature. It can 
be seen in Fig. |3that the oscillations are significant for 
UbT lower than a threshold of the order of Ae. For higher 
temperatures, the thermal fluctuations erase the struc- 
ture of superposed peaks in Q, what leads in turn to 
the disappearance of the oscillations in the a exponent. 
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FIG. 3: Log-log plot of the conductance as a function of the 
temperature. 

This agrees with the experimental observations, which 
show the existence of an inflection point T* in the log- 
log plots of conductance versus temperature, for values of 
Vg corresponding to peaks in the plot of a. As reported 
in Ref. 8, the conductance keeps following a power-law 
behavior above the temperature T*, but with an expo- 
nent significantly reduced with respect to the value at 
the peak. 

We have checked that, for values of £_f corresponding 
to peaks in the plot of a, an inflection point actually 
exists in the log-log plot of the conductance. This is 
obtained in our framework from the TDOS, by trading 
the energy variable in n{e) for the corresponding thermal 
energy ksT. The resulting dependence on temperature 
has been represented in Fig. 13 The value of the crossover 
temperature corresponds to a thermal energy of the order 
of Ae, which is consistent with the value of T* « 30 K 
that has been found in the experiment. 

We have thus achieved a description that accounts for 
the experimental features in the zero-bias conductance 
reported in Ref. 8. The discussion carried out in that 
paper highlighted the difficulty of finding a theoretical 
framework encompassing both the a oscillations and the 



existence of the crossover temperature in the conduc- 
tance. We have shown, however, that it is possible to 
extend the Luttinger liquid description of ID electron 
systems by incorporating the effects of the discrete single- 
particle spectrum, thus producing a picture that is rele- 
vant to the mentioned experimental instance. We have 
seen that the actual size of the period for the a oscil- 
lations can be explained in our approach, in terms of 
the level spacing in the quantized spectrum. The es- 
timated amplitude of the oscillations is also consistent 
with the values measured experimentally. Our theoret- 
ical construction accounts quite naturally for the exis- 
tence of a crossover in the zero-bias conductance, keeping 
track of its power-law behavior both below and above the 
crossover temperature. 

To conclude, we have developed a picture that ap- 
plies to an intermediate regime between that of Coulomb 
blockade and the Luttinger liquid behavior. In the case 
of a MWNT with a length of a few microns, for instance, 
the finite size of the system leads to a spacing between 
single-particle levels which is of order ~ 0.1 meV. Thus, 
for temperatures which are ^ 1 K we have the Coulomb 
blockade regime, with the characteristic oscillations in 
the conductance as the gate voltage is varied. For tem- 
peratures well above that value, the conductance shows a 
power-law suppression as a function of the temperature, 
typical of Luttinger liquid behavior. We have shown that 
the intermediate regime can be characterized by a power- 
law behavior of the conductance, but with an exponent 
that oscillates as a function of the gate voltage, in agree- 
ment with observations reported in Ref. This reassures 
the robustness of the Luttinger liquid picture for the de- 
scription of the CNs, down to the low-energy scales where 
the collective excitations of the electron system begin to 
fade away. 
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